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Abstract 



A full analytic analysis of the effects of CP violating phases on the event rates in the 
. direct detection of dark matter in the scattering of neutralinos from nuclear targets 



is given. The analysis includes CP violating phases in softly broken supersymmetry 
in the framework of the minimal supersymmetric standard model (MSSM) when 
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generational mixings are ignored. A numerical analysis shows that large CP violating 
phases including the constraints from the experimental limits on the neutron and the 



electron electric dipole moment (EDM) can produce substantial effects on the event 

rates in dark matter detectors. 
Supersymmetric theories with R parity conservation imply the existence of a lowest mass 
supersymmetric particle (LSP) which is absolutely stableR In supergravity mo dels! with 
R parity invariance one finds that over a majority of the parameter space the LSP is the 
lightest neutralino, and thus the neutralino becomes a candidate for cold dark matter (CDM) 
in the universe. A great deal of dark matter exists in the halo of our galaxy and estimates 
of the density of dark matter in our solar neighborhood give densities of 0.3GeV/cm 3 and 
particle velocities of ~ 320kms~ 1 . One of the interesting suggestions regarding the detection 
of dark matter is that of direct detection via scattering of neutralinos from target nuclei in 
terrestrial detectorsi. Considerable progress has been made in both technology of detection^ 
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as well as in accurate theoretical predictions of the expected event rates in detectors such 
as Ge, Nal, CaF 2 , and Xe@~i. 

In this paper we discuss the effects of CP violating phases in supersymmetric theories on 
event rates in the scattering of neutralinos off nuclei in terrestrial detectors. Such effects are 
negligible if the CP violating phases are small. Indeed the stringent experimental constraints 
on the EDM of the neutron^ and of the electronic would seem to require either small CP 
violating phases^ or a heavy supersymmetric particle spectrum@, in the range of several 
TeV, to satisfy the experimental limits on the EDMs. However, a heavy sparticle spectrum 
also constitutes fine tuning!! and further a heavy spectrum in the range of several TeV will 
put the supersymmetric particles beyond the reach of even the LHC. Recently a new possi- 
bility was proposed&il!, i.e., that of internal cancellations among the various components 
contributing to the EDMs which allows for the existence of large CP violating phases, with a 
SUSY spectrum which is not excessively heavy and is thus accessible at colliders in the near 
future. CP violating phases 0(1) are attractive because they circumvent the naturalness 
problem associated with small phases or a heavy SUSY spectrum. The EDM analysis of 
Ref.0 was for the minimal supergravity model which has only two CP violating phases. The 
analysis of Ref.0 was extended in Ref.0 to take account of all allowed CP violating phases in 
the Minimal Supersymmetric Standard Model (MSSM) with no generational mixing. This 
extension gives the possibility of the cancellation mechanism to occur over a much larger 
region of the parameter space allowing for large CP violating phases over this region. Indeed 
a general numerical analysis bears this out0. 

Large CP violating phases can affect significantly dark matter analyses and other phe- 
nomena at low energ y&M A detailed analysis in Ref.E3 shows that large CP violating 
phases consistent with the cosmological relic density constraints and the EDM constraints 
using the cancellation mechanism are indeed possible. CP violating phases affect event rates 
in dark matter detectors and a partial analysis of these effects with two CP phases and with- 
out the EDM constraint was given in RefM Thus, currently there are no analyses where 
the effect of large CP violating phases on event rates and the simultaneous satisfaction of 



the EDM constraints via the cancellation mechansim are discussed. Further, the previous 
analyses are all limited to two CP violating phases while supergravity models with non- 
universalities and MSSM can possess many more phases^!. In this paper we give the first 
complete analytic analysis of the effects of CP violation on event rates with the inclusion 
of all CP violating phases allowed in the minimal supersymmetric standard model (MSSM) 
when intergenerational mixings are ignored. We then give a numerical analysis of the CP 
violating effects on event rates with the inclusion of the EDM constraints. It is shown that 
while the effects of CP violating phases on event rates are significant with the inclusion of 
the EDM constraints, they are not enormous as for the case when the EDM constraints are 
ignored. 

We discuss now the details of the analysis in MSSM. The nature of the LSP at the 

electro- weak scale is determined by the neutralino mass matrix which in the basis (B, W, 

Hi, H 2 ), where B and W are the U(l) and the neutral SU(2) gauginos, is given by 
/ |mi|e l€l — M z sin 8 W cos (3e~ %xi M z sin 9 W sin (3e~ 

\m 2 \e^ 2 M z cos 6 W cos fie~ lXl — M z cos 6 W sin fie 

-M z sin6 w cosfie- ixi M z cos 6 W cos fie~ iX2 -|/i|e^ 

V M z sin 6 W sin fie~ ixi -M z cos 6 W sin fie~ iX2 ~M e< *" 

(1) 

Here M z is the Z boson mass, 9w is the weak angle, tan/3=|t> 2 /fi| where Vi =< Hi > 
= \vi\ e %Xi (i=l,2) where H 2 is the Higgs that gives mass to the up quarks and H\ is the 
Higgs that gives mass to the down quarks and the leptons, \x is the Higgs mixing parameter 
(i.e. it appears in the superpotential as the term fiHiH 2 ), rh\ and m 2 are the masses of 
the U(l) and SU(2) gauginos at the electro-weak scale with ^ and £ 2 being their phases. 
The neutralino mass matrix of Eq.(l) contains several phases. However, it can be shownlil 
that the eigenvalues and the eigenvectors of the neutralino mass matrix depend on only 
two combinations: 6= +X2+d^ and A£=(£i — £2)- The neutralino matrix can be 

diagonalized by a unitary matrix X such that 

X T M x oX = diag(m x o, m x o, m x o, m x o) (2) 



(4) 



We shall denote the LSP by the index so that 

X° = X* B + X* 20 W + X^H, + X* 40 H 2 (3) 

The basic interactions that enter in the scattering of the LSP from nuclei are the neutralino- 
squark-quark interactions in the s channel and the neutralino-neutralino-Z(Higgs) interac- 
tions in the cross channel. The squark mass matrix M| involves both the phases of ji and 

of the trilinear couplings as given below 
' M? + m q 2 + M 2 z {\-Q q sin 2 6 W ) cos 2(3 m q (A* q m - fiR q ) 

m q (A q m - /I* R q ) M? + m 2 + M 2 Q q sin 2 9 W cos 2(3 

Here Q u = 2/3(— 1/3) for q=u(d) and R q = v-i/v^Vz/v*) for q=u(d), and m q is the quark 
mass. The squark matrix is diagonalized by D qi j such that 

D\M\D q = diag(M| 1 , M| 2 ) (5) 

As mentioned in the introduction the relative velocities of the LSP hitting the targets 
are small, and consequently we can approximate the effective interaction governing the 
neutralino-quark scattering by an effective four-fermi interaction. We give now the result of 
our analysis including all the relevant CP violating effects in a softly broken MSSM. The 
effective four fermi interaction is given by 

£e// = X1^15XQY(APl + BP R )q + Cxxm q qq + Dxif>X™>qTihq 

+ExHbX'm q qq + Fxx^q^l^ (6) 

where our choice of the metric is i]^ u = (1, —1, —1, —1). The deduction of Eq.(6) starting 

from the microscopic supergravity Lagrangian is given in Appendix A. Here we give the 

results. The first two terms (A, B) are spin-dependent interactions and arise from the Z 

boson and the sfermion exchanges. For these our analysis gives 

r ,2 \r< |2 \n' |2 



B = -l4 l|A '»' 2 " W Wft, + 4 (4-%) + 4(M|-'%) <8) 



where 



C q L = V2{a q0 D q ii - 7gO-Dg2li 



C q R — V2([3 q oD q u — 5 q oD q 2l) 



C' qL = V2{a q0 D ql2 - 7 g o-D: 



q22) 



(9) 
(10) 

(11) 



C' qR = V2(P q0 D ql2 - 5 q0 D q22 ) (12) 
and where a, /3, 7, and 6 are given byil 

a u(d)j - ~ : — w, [i-o) 

2mwsm/3(cos/3) 

Pu(d)j — eQu(d)jX[*j + — ^ — X' 2 *(T 3u ^ — Q u (d) sin 2 6 W ) (14) 

- *n y> 9 ® u ^ sin2 9w v cik\ 

lu(d)j - eQ u{d)j X 1:i — — X 2j (15) 

, -9m u{d) X* (3):j 

u(d)j ~ 2m w sm(3(cosP) [ ' 
Here g is the SU (2) L gauge coupling and 

X[j = X X j cos9 w + X 2j sin8 w (17) 



X' 2j = —Xij sin8 w + X 2j cos9 w (18) 

The effect of the CP violating phases enter via the neutralino eigen-vector components X^ 
and via the matrix D qi j that diagonalizes the squark mass matrix. 

The C term in Eq.(6) represents the scalar interaction which gives rise to coherent scat- 
tering.. It receives contributions from the sfermion exchange, from the CP even light Higgs 
(h°) exchange, and from the CP even heavy Higgs (H°) exchange. Thus 
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C = C ' j + Cyfi + CffO 

where 

^ , , , g 2 cos afsma) 

y ' ; v '4M W M* sin /3 (cos/3) 

g 2 sin a (cosa) 
= AM W M 2 H0 sin (3(cos(3) Rep 

Here ct is the Higgs mixing angle, (u,d) indicate the flavor of the quark involved 
scattering, and a and p are defined by 

a = X^X^q — tan#vi^io) cos a + X^X^q — tan 9wX* ) sin a 

p = — X^X^q — tam9wX* ) sin a + X^X^q — tan 9wX^ ) cos a 
The last three terms D,E and F in eq.(6) are given by 



D(u,d) = CAu.d) H k Reuj 

/v AM W m\ 



q 2 , , .cosalsina) Ima sina(cosa) Imp^ 

E(u, d) = TAu, d) + -f — - +— Jr^r- + —77 7F7 — 

/v ' 4M w l v ' sm/3 (cos/3) m 2 sm /3 (cos /3) m 2 H o 



, ,. //- cotf3(tanf3) 

Flu, d) = T-Au, d) + —— ^ Imuj 

fK AM W m 2 A0 



where A is the CP odd Higgs and where u is defined by 

u = -X 4 * (X 2 * - taa6 w X* ) cos/? + X 3 * (X 2 * - t&n9 w Xl ) sin (3 

and 
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In the limit of vanishing CP violating phases our results of A, B and C limit to the result 



Ref.H in their equations (2a, 2b, Al). Further, in the same limit of vanishing CP violating 
phases our results go to those of Ref.i with an overall minus sign difference in the three 
exchange terms, i.e., the Z, the sfermions and the higgs terms, (see Appendix B for details). 



limited to the case of two CP violating phases and it gave the analytic results for only 
some of the co-efficients in the low energy expansion of Eq.(6). The terms E and F given 
by Eqs.(26) and (27) are new and vanish in the limit when CP is conserved. The term D 
given by Eq.(25) is non- vanishing in the limit when CP phases vanish. However, this term 
is mostly ignored in the literature as its contribution is suppressed because of the small 
velocity of the relic neutralinos. In fact the contributions of D,E and F are expected to be 
relatively small and we ignore them in our numerical analysis here. 

For the computation of the event rates from nuclear targets in the direct detection of 
dark matter we follow the techniques discussed in Ref.i and we refer the reader to it for 
details. We give now the numerical estimates of the CP violating effects on event rates. First 
we consider the case when the EDM constraint is not imposed. In Fig.l we exhibit the ratio 
R/R(0) for Ge where R is the event rate with CP violation arising from a non-vanishing 9^ 
and R(0) is the event rate in the absence of CP violation. The figure illustrates that the 
effect of the CP violating phase can be very large. In fact, as can be seen from Fig.l the 
variations can be as large as 2-3 orders of magnitude. However, as noted above the EDM 
constraint was not imposed here. Next we give the analysis with inclusion of the EDM 
constraints. For this purpose we work in the parameter space where the cosmological relic 
density constraint and the EDM constraints are simultaneously satisfied and we compute the 
ratio R/R(0) for Ge in this region. Specifically the satisfaction of the relic density and the 
EDM constraints is achieved by varying the parameters of the theory. The satisfaction of the 
EDM constraint is achieved through the cancellation mechanism discussed in Refs. 0>0. The 
result of this analysis is exhibited in Fig.2. Here we find the range of variation of R/R(0) 




To compare our 



results to those of referenced we note that the analysis of 




was 
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with dfj, to be much smaller although still substantial. Thus from Fig. 2 we find that the 
variation of R/R(0) has a range of up to a factor of 2 over most of the allowed parameter 
space satisfying the relic density and the EDM constraints in the regions of the parameter 
investigated. This variation is substantially smaller than the one observed in Fig.l when the 
EDM constraints were not imposed. 

In conclusion, we have given in this paper the first complete analytic analysis of the effects 
of CP violating phases on the event rates in the direct detection of dark matter within the 
framework of MSSM with no generational mixing. We find that the CP violating effects 
can generate variations in the event rates up to 2-3 orders of magnitude. However, with the 
inclusion of the EDM constraints the effects are much reduced although still significant in 
that the variations could be up to a factor ~ 2 as seen from the analysis over the region of 
the parameter space investigated. Of course the parameter space of MSSM is quite large 
and there may exist other regions of the parameter space of MSSM where the CP violating 
effects on event rates consistent with the EDM constraints are even larger. 
This research was supported in part by NSF grant PHY-96020274. 
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FIGURES 
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Fig. 1. Plot of R/R(0) for Ge as a function of the CP violating phase 6^ in the MSSM case when 
tan/3=2, mo=100 GeV, |A)|=1 f° r the cases when the gluino mass is 500 GeV (solid curve), 600 
GeV (dotted curve), and 700 GeV (dashed curve) 
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Fig. 2. Plot of R/R(0) for Ge as a function of the CP violating phase 6 ^ for the parameter space 
discussed in the text. The parameter space spans regions satisfying the relic density and the EDM 
constraints obtained by varying other parameters in the theory. 
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1. Appendix A 



In this appendix we give a derivation of the four fermi neutralino-quark effective Lagrangian 
with CP violating phases given in the text. 

A. squark exchange terms 

From the microscopic Lagrangian of quark-squark-neutralincB 

-C = q[C qL P L + C qR P R ] X <h + q[C' qL P L + C' qR P R }xq~2 + H.c. (30) 
the effective lagrangian for q — x scattering via the exchange of squarks is given by 
c *ff = M 2 z M 2 x[c; l Pr + C; R P L }qq[C qL P L + C qR P R ] X 

91 X 

+ M2 _ M 2 X[C;'lPr + C; R P L ]qq[C' qL P L + C qR P R } X (31) 
ivi q2 ivi x 

We use Fierz reordering to write the Lagrangian in terms of the combinations xxqq-i 
XlbXqisq, X7 m 75X?7m?>X7''75X?VM j XTtsXW and xxqiM- For this purpose, we define the 
16 matrices 

r A = {l,7 ,^,2 7 75 ,y 75 , 75 , W °\^} : i,j = 1-3 (32) 
with the normalization 

tr(T A T B ) = 46 AB (33) 
The Fierz rearrangement formula with the above definitions and normalizations is 
(u 1 r^ 2 )(u 3 r B n 4 ) = £ F A E(uiT C u 4 )(u 3 T D u 2 ) (34) 

C,D 

where Uj are Dirac or Majorana spinors and 

FcE = -(+)^tr(T c T A T D T B ) (35) 

and where the +ve sign is for commuting u spinors and the -ve sign is for the anticommuting u 
fields. In our case we have to use the -ve sign since we are dealing with quantum Majorana 
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and Dirac fields in the Lagrangian. We give below the Fierz rearrangement of the four 
combinations that appear in Eq.(31) above: 

11 1 u 

xqqx = -^xxqq - -^xi^xQim + txttsx^vm 

i_ i_ 1 _ 

X15QQX = -^CTlhXTi^ - ^xxwsq - ^xmxQQ 

1 „ 1 1 

XQQ15X = —^XTlhXTIiiq. ~ -^XXQlbQ ~ ^XlbXQQ 

1 _ 1 _ _ 1 _ 

XimisX = -jXXQQ - 4X75X^754 - 4X7^75X^75? (36) 

In the above we have used the metric rj^ = (1,-1,-1,-1), and we also used the fact that 
X is a Majorana field so that X7^X = an d X a nuX = 0- By rearranging the terms to be in 
the form of C e ff of Eq.(6) we can find out directly the squark contributions to A, B, C, D, E 
and F terms as given in the text of the paper. 



B. Z boson exchange 

From the q — Z — q and x ~ % — x interactions in Eqs. (c62) and (c87a) of Ref.@, we obtain 
the following effective Lagrangian for the q — x scattering via Z-exchange: 

C *ff = AM2 g2 2« [l^3o| 2 - \X A ,\ 2 ]qY[d qL P L + d qR P R ]qxi»lsX (37) 

4iW^ COS u\v 

where d q i = T 3 l — e q sin 2 9w and d g n = —e q sin 2 6w- From Eq.(37) we can read off the 
contribution to A and B from the Z exchange. These contributions are given in Eqs. (7) and 



C. Higgs exchange terms 

Higgs exchange will contribute to C,D,E and F terms. From the interaction Lagrangian 
of Chxx an d ^Hqq m Eqs. (4.47) and (4.10) respectively of Ref.0, one can get the effective 
Lagrangian for q — x scattering via h°, H° and A exchanges. In our formalism we use h°, 
H° and A for the light, heavy and CP-odd neutral higgs. There are six contributions: three 
higgs exchange terms for the up flavor and three for the down flavor. To illustrate we choose 
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the up quark scattering with \ y i a the exchange of the heavy higgs H° (H± in the notation 
of RefH): 



where 



m m 



J ho = -§ cosaxiQloPL + Qm P R)X 



J m = | sinax(^oo p i + s oo p r)x 
gm u sin a _ 

^ = ~2M^P UU (39) 
where Qqq, Sqq are as defined in Ref.il. Defining p by 

P = Qoo cos a ~ Sqo s i n a (40) 

we get the H° contribution to £> e ff 

g 2 m u sin a , , _ 
^ = 4M w m^ sin/ e(p) ^ 
, ig 2 m u sina 

—-Im{p)xi^xuu (41) 



4M w m 2 HQ sin /3 

From Eq.(41) we can read off directly the contributions Cho(u) and Eho(u) as given by 
Eqs.(22) and (26). 
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2. Appendix B 

Here we compare our results with those of Ref.i which is in the limit of no CP violation, 
no sfermion mixing and no heavy Higgs. In the limit of no CP violation and no sfermion 
mixing Cl,C' l , C r , C' r given by Eqs.(9-12) in the text reduce to the following: 

C L = V2a q0 , C' L = - V2~/ q0 

C R = V2(3 q0 , C' R =-V25 q0 (42) 

To express the above in the notation of Ref.i we set X± = j3, X^ = a, X^ = 5, and 
Xl = 7 keeping in mind that the Hi and H 2 of Ref.i are defined oppositely to our notation. 
Using the above along with Eqs. (13-16) we find for T 3 = |(— ~) 

2 _ ™ 2 u(™ 2 dh 2 ($ 2 ) 



\Cr 



» 2 M) 

\C' L \ 2 = 2{ 9l l -Y R (5f 



\C R \ 2 = 2(ag 2 T 3 + Pg 



Y, 



L \2 



2 



\ n , 12 _ ™l(™ 2 dh 2 (S 2 ) (a o\ 

1 * l ~ WA) { } 

where Yl is the hypercharge, V\ =< Hi >, and v 2 =< H 2 > and where Hi and H 2 are in 
the notation of Ref.i. Further using the identity 

9 2 1 
2— (T 3 - e q sin 2 9 w ) = -{g x sin6 w + g 2 cos9 w )(-Y L gisin9 w - T 3 g 2 cos9 w ) (44) 

cos 2 9 w 2 

where the left hand side of Eq.(44) is written in the form used in Ref.i, we can express A 
and B in the limit of no CP violation and no sfermion mixing as follows: for T 3 = ~(— |) 

(7 2 - <5 2 ) 1 
A = — (gism9 w + g 2 cos9 w )(-Y L giszn9 w - T 3 g 2 cos9 w ) 

z 

^92T 3 + (3gi^f K(m 2 d h 2 (S 2 ) 



Ml) ±{Mh-Ml)v 2 M) 

(l 2 ~ S 2 ) 1 
B = 2 — {gisinOw + g2Cos9 w )(-Y L gism9 w - T 3 g 2 cos9 w ) 

~ z ~~ 



+ 2(M| R - M*) + 4(Ml - M^v\(yl) 1 ° j 
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To compare our C term with that of RefH we again go to the limit of vanishing CP violating 
phases, assume no sfermion mixing, and in addition ignore the heavy Higgs exchange con- 
tribution (i.e., the term C#o of Eq.(22) in the text). Then using similar notational changes 
as above we find that our C under the approximations made in RefJ is given by 

_ g\ (—cosa)(sina), _ n . . . 

AM w m 2 h0 (sin/3)(cosf3) ' 
92 ,ag 2 T 3 + p gi ^- (3 9l ^ N 7 ($) m 1,1 



AM W K M\ L -Ml M\ R - sinp(cospy Ts 2^ 2^ ^ 

Comparing our results for A,B and C with those of Ref.i we find that our Z, sfermion and 
Higgs exchange terms have an overall minus sign relative to those of Ref.ll. 
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